
######################################################################
######################################################################
######################################################################
######  This file contains code to run all r-scripts            ######
######      to replicate all results in                          #####
######     in Betz, Cook, Hollenbach 2019                       ######
#### Bias due to network misspecification under spatial dependence ###
######################################################################
######################################################################

#### load all necessary packages
library(tidyverse)
library(ggplot2)
library(here)
library(xtable)
library(purrr)
library(sp)
library(spdep)
library(MASS)
library(spatialreg)
library(foreach)
library(doParallel)
library(haven)
library(tictoc)

sink(here("..", "results", "log_sessionInfo.txt"))

sessionInfo()
sink()

#### applied example mentioned in text
source(here("examples.R"), echo = T)

#######
###### running simulations, output saved in data folder
source(here("sim_sar_minmax.R"), echo = T)
source(here("sim_sar_row.R"), echo = T)
source(here("sim_slx_minmax.R"), echo = T)
source(here("sim_slx_row.R"), echo = T)
#####
#### now create plots and tables, output saved in results folder
source(here("SAR_plots.R"), echo = T)
source(here("SLX_plots.R"), echo = T)
